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Analytical solution of one dimensional hydrodynamical model is derived, where phase transition 
from the QGP state to the hadronic state is effectively taken into account. The single particle 
rapidity distribution of charged n mesons observed in relativistic heavy ion collisions is analyzed by 
the model. Space-time development of the fluid is also investigated. 
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I. INTRODUCTION 

It is widely recognized that at the initial stage of rel- 
ativistic heavy ion collisions, the quark-gluon plasma 
(QGP) is formed, and it breaks up to thousand of 
hadrons in the final states. It is also considered that hy- 
drodynamical approach is one of the effective tools [ll-Q 
to analyze such multiple particle production processes. 

In most of numerical calculations of three dimensional 
hydrodynamical models, the phase transition from the 
QGP state to the hadronic state is taken into account 0- 
6j. For reviews of the recent works, see Q- 

Recently, exact solutions of one dimensional hydro- 
dynamical model are investigated [§-[Tll|. However, the 
phase transition from the QGP state is not introduced 
to the analytical formulations. 

The hydrodynamical model is consisted of the energy 
momentum conservation, the baryon number conserva- 
tion (or entropy conservation) and the equation of state. 
The phase transition from the QGP state to the hadronic 
state can be expressed by the change of the velocity of 
sound from that of the perfect fluid, cq = l/v3, to some 
constant value, c s , at critical temperature. We assume 
that shear viscosity, bulk viscosity and heat conductiv- 
ity are negligibly small in the hadronic state. We also 
assume baryon number can be neglected. 

One dimensional hydrodynamical model of the perfect 
fluid can be solved with an appropriate initial condition. 
If the change of the velocity of sound is taken into the 
model, we can formulate the hydrodynamical model in- 
cluding the phase transition from the QGP state to the 
hadronic state. 

In section [TTl one dimensional hydrodynamical model 
following to Landau's approach is reviewed, and a solu- 
tion for potential x is derived under an idealized initial 
condition. In section IIII1 one dimensional hydrodynam- 
ical model including the phase transition from the QGP 
state to the hadronic state is formulated. In section HVl 
the single particle distribution under the Cooper- Frye ap- 
proach is shown. The single particle rapidity distribution 
observed in relativistic nucleus-nucleus (AA) collisions is 
analyzed in section [V] The space-time evolution of fluid 
element is also investigated there. The final section is 
devoted to concluding remarks. 



II. ONE DIMENSIONAL HYDRODYNAMICAL 
MODEL 

A. Equations for hydrodynamical model 

The hydrodynamical model proposed by Landau 
is composed from the energy-momentum conservation of 
the perfect fluid and the equation of state. 

The energy-momentum conservation is given by, 

JLt^ = o, ^ = 0,1, (1) 

where T^„ denotes the energy momentum tenser of the 
perfect fluid, 



Tpv = (e + p)u^u u - pg^ y 



(2) 



In Equation , e denotes the energy density and p the 
pressure of a fluid element. The velocity of the fluid 
element is denoted by u M , which satisfies, u^u^ — 1, and 
3^ = diag(l,-l). 

The equation of state is expressed as a function of the 
energy density s to the pressure p. It would be different 
whether the fluid is in the QGP state or in the hadronic 
state. At present we assume that the relation, dp/de = 
cq 2 , holds, where cq is a positive constant. 



In addition, thermodynamical relations, 



- P 



Ts, 



Tds, 



(3) 



are used, where T denotes the temperature and s is the 
entropy density of the fluid element. Equation © holds 
whether fluid is in the QGP state or in the hadronic state. 

Projection of Eq.([T]) to the direction of u M gives the 
entropy conservation, 



d{su v ) 
dx v 



= 0. 



(4) 



After Eq.([T]) is projected to the direction perpendicular 
d(Tu v 



to Up, we obtain 



dT 

dx^ dx^ ' ^ ' 

Rapidity r\ of the fluid element is defined by, 

(u ,^ 1 ) = (cosh 77, sinhr;). 

Then, both equations in Eq.([5]) reduce to 

d (3 
— (Tsinhr/) + — (T cosh 77) = 0. 
at ox 
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(5) 



(6) 



(7) 
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From Eq.([7]), there is a function, (j>, which satisfies, 

^ = Tcosh77, ^ = -Tsinhrj. (8) 
at ox 

By the use of the Legendre transform, 

d\ = d{4> — tT cosh 77 + xT sinh 77) , 
we obtain the equations for the potential \ [U, H| ; 

dx 

dT 
ldx 
Tdri 



-t cosh rj + x sinh 77, 
-t sinh 77 + x cosh 77. 



(9) 



Then, space-time variables of the fluid element, t and 
x, are given respectively asC 

1 / dy dx \ 

t = — e" ( 77— cosh 77 + — sinh 77 ) , 



T \duj 



Or] 



x = -U" sinh cosh ry), (10) 

where To is the initial temperature of the fluid and uj = 
ln(T /T). 

The entropy conservation, Eq.QCis rewritten as, 
d d 

— (scosh?/) + — (ssinhr/) = 0. (11) 

After changing the variables from x, rj to 77, T in Eq. (fTT|) . 
we obtain, 

T ds /d 2 x r dx\ _ T 2^X = n 
sdTKdr] 2 dT) dT 2 



B. Solution for potential \ 

Using the relation (T / s)ds / dT = cq 2 and the variable, 
u> = — ln(T/T)), we obtain the partial differential equa- 
tion for the potential X) 



d 2 x 2 ®X ^ d 2 x g 
duj 2 ° duo co 2 drf ' 

B - 1 " c ° 2 -l 



(12) 



Equation (fl"2)) is called the equation of telegraphy. 

The initial condition for Ea. (lT2l is taken with a con- 
stant Qo as, 



dx 

duj 



j=0 



G(v) = QoS(v)- (13) 



Hydrodynamical model is applicable mainly to the cen- 
tral region [l2j|. As can be seen from Eq. (fT7|) . the main 
contribution to the fragmentation comes from the func- 
tion g(rf). Therefore, we assume that g{rj) = 0. At the 
initial stage, fluid would be formed in the very small re- 
gion of rapidity space, and we put G(rf) is proportional 
to the delta function in the rapidity space as an idealized 
case 1131. 



Introducing the new variable Xi by, 
X (v,u,)=xi(v^y° u , 
we have the partial differential equation for xi as ; 



d 2 Xi a 2 



Po Xi 



1 d 2 X i 

ou z co z or} 2 - 

and the initial condition for Eq. (fT5")) as 

Xi|o;=o = g(v) = 0, 
dxi 



= 0, 



(14) 



(15) 



did 



= G(r 1 )~B g(r ] ) = Q S( V ). 



(16) 



The solution for Xi is given in general as [14] . 

Xifow) = ^{g(v + ^/c )+g(r]-u/c )} 



£0 
' 2 



j/co 



dz 



—oj/co 



Bow r" /co 



[G(z + ri)-Pog{v)} 



' LU 2 — Cq 2 Z 2 



4/ CO 



dz g(z + 77) 



h{Mu 2 - tfz 2 ) 



\Jlo 2 — co 2 z 2 



It reduces to 



Qoco 



I (Bq \Juj 2 - CQ 2 !] 2 



(17) 



(18) 



for |r7| < cj/c . 

Then the solution for Eq. (TT2")) under the initial condi- 
tion (p~3]) is given by 



x(v,u) 



co 2 ?? 2 



(19) 



III. HYDRODYNAMICAL MODEL INCLUDING 
PHASE TRANSITION 

A. Introduction of phase transition 

In our formulation, the phase transition from the QGP 
state to the hadronic state is expressed by the change of 
the velocity of sound. The evolution of fluid is described 
by the parameter u> = ln(To/T), where To is the initial 
temperature. The fluid formed just after the collision of 
nuclei, is assumed to be in the QGP state. The velocity 
of sound in it is denoted by Co (cq = 1/ Vo) . It expands 
from co = to uj — u) c — 0, where w c = ln(To/T c ) > 
with the critical temperature T c . The phase transition 
from the QGP to the hadronic state occurs at u> = ui c . 
From ui = uj c + to to = ujf, where uif = ln(To/T/) > ui c 
with the freeze-out temperature Tj, the fluid is in the 
hadronic state. The velocity of sound in it is assumed to 
be constant and denoted by c s (0 < c s < cq). 
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The equation of state is assumed to be, 

c 2 £ - (1 + co 2 )B, for < oj < u) c , 
for oj c < oj < ujf, 



P 



c s e, 



(20) 



where B is the bag constant of hadrons. From Eq. (f20)) . 
the velocity of sound is given by 



Co, for < oj < oj c , 
c s , for ou c < oj < oj f . 



(21) 



As is shown in Figure [TJ the velocity of sound is discon- 
tinuous at oj = oj c . 
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FIG. 1. the velocity of sound as a function of oj from oj — 

to OJ = OJf 



B. Potential x(v,^) in the QGP fluid 

We assume that the fluid in the QGP phase, where the 
velocity of sound is c = 1/ continues from oj = to 
oj = oj c — 0. Then, the partial differential equation for 
the potential x(Vi oj) is given by Eq. (fT2|k 



d 2 \ 2 ^ <9 2 X q 

doj 2 doj cq 2 drf ' 



(22) 



where /3o = (1 — co 2 )/2co 2 = 1. The initial condition for 
Eq.flm is taken as Eq.^, 

xl«=o = 9(v) = 0, ^ = G( V ) = QoS( v ). (23) 
ooj lj=o 



The solution for Eq. ([22|) under the initial condition 
23} is given by, 



QqCo 



o(/W 



oj 2 — Co 2 rj 2 



(24) 



For oj c < oj <ojf, the partial differential equation for 
the potential x is given by, 



Q 9 " Q 9 O 9 ' r 1 

ooj 1 ooj c s z or\ A 



2c, 



> 1.(25) 



For \r]\ < (oj — oj c ) /c s + a; c /coCthe solution for Eq. ([2~5)) is 
given by 



Xfau) = A(oj)I (X), 



where 



(26) 



(27) 



A = PCsy/vlnax ~ V 2 , 
Vmax = (W - OJ c )/c s + 0J c /cq. 



If c s = Co (ft = /?o)Cthe potential xiv^) in Eq.(|26 
coincides with Eq. (|24|) . which is the solution for Eq.([22 
under the initial condition ([2"3")l D 



D. Effect of phase transition on potential \(rj,oj) 

In our formulation, the velocity of sound is changed 
from Co to c s according to the phase transition from 
the QGP state to the hadronic state. This effect would 
appear in the potential x(Vi w ) m the neighborhood of 

oj = 0J C . 

At oj = ou c — 0, where the fluid is in the QGP state, the 
potential is given from Eq. ([24|) as, 



X (r), oj c -0) = ^V^/o (po V"c 2 - c V ) • (28) 

At oj = oj c + 0, where the fluid is in the hadronic state, 
it is given from Eqs. (l2l)l) and (|2"T1) as, 

X(v,u>c + Q) = ^V^/o( — V^ 2 - C qV)-(29) 
2 V co / 

From Eqs.([2"g]) and ([Ml), we obtain, 

x(^-0) /o(/3oV^c 2 -coV) 

The potential x(r;, oj) becomes discontinuous at oj = oj c 
due to the phase transition. This fact will affect the 
space-time behavior of fluid element. 



C. Potential x(r},oj) in the hadronic fluid 



IV. SINGLE PARTICLE RAPIDITY 
DISTRIBUTION 



In the hadronic phase, the fluid expands from oj = 
oj c + to oj = ojf . During the expansion of the hadronic 
fluid, with the velocity of sound is c s . 



In order to analyze the rapidity distribution of ob- 
served particles, the Cooper-Frye approach [TEj is used. 
Then, the single particle rapidity distribution dn/dy of 



4 



particles with mass m and transverse momentum px in nucleus-nucleus collisions is given byC 

dn ttRa 2 f ( dx dt 

^ = ^JS coshy d^- smhy T v 

rriT 



r ut xTT^n — rdr}d p T , (31) 
exp[my cosn(y — T))/ 1 /J — 1 

where Ra denotes the radius of colliding nuclei and 
tut = V 'pt 2 + m 2 denotes the transverse mass of ob- 
served particles. Space-time variables, t and x, in Eq. (|31[) 
is given by Eq. fTU)) . Under the assumption that the 
Bose-Einstein distribution can be approximated by the 
Maxwell-Boltzmann distribution, Eq. (j3"Tj) is written as, 



dn RA 2 Tf 2 



4tt 
/ m 2 



dy 



dt\ 
2m 



d 



T f 2 T f cosh(y - 77) cosh 2 (y - n) 



tanh(y-7j) + c s 2 ^-(x + ^) 



-m cosh(y— rj) /Tt 



(32) 



If the term multiplied by tanh(y — rj) is neglected in the 
square bracket on the right hand side of Eq. ([32l) . the 
single particle rapidity distribution coincides with the 
Milekhin's approach [16j,ll7| except for the normalization 



factor. 

By the use of Eqs. (|2"oT) and (|2"7|) . the single particle 
distribution is given by, 



dn Ra'T 2 



dy 



4tt 



WcfA^,) f~ A,[^{5S£|^/ (A) + (^±l-25S*^)/ 1 (A)}ta„h( !/ -,) 



-{( 



1 + /3 (/3c s ?7 mQX ) 2 \ {) (n. 



A 2 
2m 



A V c s 
2 



1 - 2 



2 



,T/ 2 Tf cosh(?/ — 77) cosh 2 (y-r/) 



A 2 

e -mcosh(y— »7)/T/ 



A 2 
i"i(A)} 



(33) 



Space-time variables of the fluid element are written from 
Eqs.||TDJ} and J2BJ) as, 



t = y a ^ K /o(A) + ^-^-fiW} coshr ? 



-(/3c s ) 2 |/ 1 (A)sinh7 7 



(34) 



^A(u)c" \p{l Q (\) + Pc s ^h(\)\ sinhr? 
J L L A J 



-(/3c s ) 2 |/ 1 (A)cosh ?7 



DATA ANALYSIS 



(35) 



The single particle rapidity distribution of charged tt 
mesons in Au-Au collisions at 200 AGeV at RHIC [3 
is analyzed by Eq. (|33l) . In the analysis, the initial tem- 
perature Tq is fixed at 0.95 GeV, the critical tempera- 



ture T c at 0.18 GeV and the freeze-out temperature Tf 
at 0.12 GeV. The nuclear radius Ra is parametrized as, 
Ra = 1.2 x A 1 / 3 , where A denotes the mass number of 
nucleus, and A = 197 for Au. Then, two parameters, 
Qo and c s , remains in our formulation. Estimated pa- 
rameters are shown in Table 1, and the result on the 
single particle rapidity distribution of charged tt mesons 
is shown in Fig. [2] 



TABLE I. Estimated values of parameters Qo and c s 2 from 
Au+ Au -> (tt+ +tt~)/2 + X 



Qo 



(4.78 ± 1.15) x 1(T 2 0.117 ± 0.005 41.4/12 

The space-time evolution of fluid element at 
w=constant in the x-t plane is shown in Fig|3] The 
profile at T = T c = 0.180 GeV is overlapped to the 
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FIG. 2. Analysis of single particle rapidity distribution of 
charged pions (tv + tt~)/2 by Eq. (|33|l . 



FIG. 4. (Color online) Space-time evolution of fluid element 
at fixed temperature from T = 0.180 GeV to T = 0.170 GeV. 



origin (x,t) = (0,0). The dotted curve denoted by 
BJ in the fig ure means the Bjorken's scaling solution, 
t = \Jt 2 — x 2 =constant [191 . 
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\Jt 2 — x 2 , in FigJS] The temperature profiles of fluid ele- 
ments for T > 0.200 GeV are abbreviated. 
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FIG. 3. (Color online) Space-time evolution of fluid element 
at fixed temperature 

In order to investigate the space-time evolution of fluid 
element around u) = uj c or T = T c , the detail of it is shown 
in FigH The profile at oj = uj c - or T = T c + is in 
the QGP state. It is calculated by Eq.([2i]). The profiles 
at uj > uj c or T < T c are in the hadronic states. Those 
are calculated by Eqs. (l2l)|) and (|2"T|) . As can be seen from 
FigHl the region between the profile at T = T c + and 
that at T = T c — corresponds to the mixed phase where 
the QGP state and the hadronic state are coexistent. 

The temperature dependence of fluid element at r\ = 
constant is shown as a function of proper time, r = 



FIG. 5. (Color online) Temperature profiles of fluid elements 
at fixed rapidity. 

The fluid element with larger rapidity in the absolute 
value is cooled faster, as can be seen from the figure. The 
fluid starts expansion from T = T at 77 = and at r = 0. 
For example, the fluid element at r\ = is cooled down 
from T = T at r = to T = T c at r = 0.238 fm/c. It 
remains at T = T c from r = 0.238 fm/c to r = 4.957 
fm/c. Therefore, the mixed phase in the neighborhood 
of 77 = is to continue about 4.7 fm/c. The fluid element 
at 77 = is cooled down to T — Tf at r = 134.9 fm/c. 

The dotted curve denoted by BJ corresponds to 
the Bjorken's scaling solution, T = T c (t c /t) Cs 
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with c s 2 = 0.117, which is normalized to T c at 
(t,T)=(t c ,T c )=(0.238,0.18) . 

VI. CONCLUDING REMARKS 

We have formulated one dimensional hydrodynamical 
model including the phase transition from the QGP state 
to the hadronic state. At first, following to the Landau's 
hydrodynamical model, the equation of the telegraphy 
for the potential xijli w ) is solved under the simplified 
initial condition in the rj — oj space with the velocity of 
sound, Co = 1/ v3- Then the solution of it in the hadronic 
state with the constant velocity of sound, c s (c s ^ Co), 
is found so as to coincide with the solution in the QGP 
state if c s — cq. 

The space-time evolution of fluid element from w = 
(or T = T Q ) to ui = ujf (or T — Tf) is calculated by our 
model. The non-zero finite region emerges between the 
profile at w = w c — (or T = T c +0) and that at w = w c +0 
(or T = T c — 0). This is due to the discontinuity of the 
potential xiVi w) at w = w c (or T — T c ) caused by the 
change of the velocity of sound from Co in the QGP state 
to c 8 (c 8 < Co) in the hadronic state at w = uj c . In 
our calculation, the fluid element at rj = freezes out at 
t = 134.9 fm/c. 

In the computer simulations of the three dimensional 
hydrodynamical models, the evolution of fluid element 
is calculated as a function of the proper time r. For 
example, the initial condition is taken as t = 0.6 fm/c 
and T = 0.36 GeV at 200 AGeV j|. 

In our one dimensional hydrodynamical model, the ini- 
tial condition is taken as r = fm/c and To = 0.95 GeV. 
The fluid element at rj = is cooled down to T = 0.36 
GeV at t — 0.036 fm/c. Therefore, the temperature of 
fluid in our calculation at t = 0.6 fm/c is not higher than 
that of the three dimensional hydrodynamical model. 



In the computer simulations at 200 AGeV, the fluid 
element expands a few hundred fm along the direction of 
colliding nuclei before freeze-out though the value of r at 
freeze-out is about 15 fm/c or so (2pj . 

The proper time of fluid element at rj = at freeze-out 
in our calculation is one order larger than the proper time 
of fluid element at freeze-out in the three dimensional 
computer simulation. In the one dimensional hydrody- 
namical model, the expansion of fluid into the transverse 
dimension is neglected contrary to the three dimensional 
calculations. Therefore, the temperature of the fluid el- 
ement especially in the neighborhood of rj — would 
decrease more slowly than that calculated by the three 
dimensional computer simulation, where the transverse 
expansion of the fluid is taken into account. However, 
the scale of expansion of fluid element in the space vari- 
able in our calculation shown in FigO is comparable to 
that along the direction of colliding nuclei in the com- 
puter simulation. 

In the recent lattice QCD calculations, it is shown that 
the crossover phase transition occurs at chemical poten- 
tial \i = [il], |2^| . Let the crossover transition take place 
in the region, (u> c — Aw/2, w c + Aw/2). In our formu- 
lation, it would be expressed by the smooth change of 
the velocity of sound from Co to c s . The sound of veloc- 
ity c v (ui) in (w c — Aw/2, w c + Aw/2) is assumed to be 
a continuous function with positive value, and to satisfy 
c„(w c — Aw/2) = Co and c„(w c + Aw/2) = c s . 

The potential x(w) in (0, w c — Aw/2) is the same with 
Eq. lfMj) . The potential x(lo) after the crossover transi- 
tion, in (w c + Aw/2, w/), would be given by the following 
equation, 

xfa,u)=A(w)Io(A), (36) 

where, 



(37) 
(38) 

I 
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A{uj) = ^pexp[/3(w - w c - Aw/2) + A,(w c - Aw/2) + 



^ = fiCsVVmax - V 2 , Vmax = (w - W c - Aw/2)/c s + (w c - Aw/2)/c + 



p v (x)dx] , /3 v (x) = — 2- 

• a -Au/2 2c v (X) 



-Au)/2 



c v (x) 



dx. 
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